Optimal approximate matrix product in terms of stable rank

We prove, using the subspace embedding guarantee in a black box way, that one can achieve the spectral norm guarantee for approximate matrix multiplication with a dimensionality-reducing map having $m = O(\tilde{r}/\varepsilon^2)$ rows. Here $\tilde{r}$ is the maximum stable rank, i.e. squared ratio of Frobenius and operator norms, of the two matrices being multiplied. This is a quantitative improvement over previous work of [MZ11, KVZ14], and is also optimal for any oblivious dimensionality-reducing map. Furthermore, due to the black box reliance on the subspace embedding property in our proofs, our theorem can be applied to a much more general class of sketching matrices than what was known before, in addition to achieving better bounds. For example, one can apply our theorem to efficient subspace embeddings such as the Subsampled Randomized Hadamard Transform or sparse subspace embeddings, or even with subspace embedding constructions that may be developed in the future. Our main theorem, via connections with spectral error matrix multiplication shown in prior work, implies quantitative improvements for approximate least squares regression and low rank approximation. Our main result has also already been applied to improve dimensionality reduction guarantees for $k$-means clustering [CEMMP14], and implies new results for nonparametric regression [YPW15]. We also separately point out that the proof of the "BSS" deterministic row-sampling result of [BSS12] can be modified to show that for any matrices $A, B$ of stable rank at most $\tilde{r}$, one can achieve the spectral norm guarantee for approximate matrix multiplication of $A^T B$ by deterministically sampling $O(\tilde{r}/\varepsilon^2)$ rows that can be found in polynomial time. The original result of [BSS12] was for rank instead of stable rank. Our observation leads to a stronger version of a main theorem of [KMST10].Comment: v3: minor edits; v2: fixed one step in proof of Theorem 9 which was wrong by a constant factor (see the new Lemma 5 and its use; final theorem unaffected

Population Dynamics on Complex Food Webs

In this work we analyse the topological and dynamical properties of a simple model of complex food webs, namely the niche model. In order to underline competition among species, we introduce "prey" and "predators" weighted overlap graphs derived from the niche model and compare synthetic food webs with real data. Doing so, we find new tests for the goodness of synthetic food web models and indicate a possible direction of improvement for existing ones. We then exploit the weighted overlap graphs to define a competition kernel for Lotka-Volterra population dynamics and find that for such a model the stability of food webs decreases with its ecological complexity.Comment: 11 Pages, 5 Figures, styles enclosed in the submissio

Shock Formation in a Multidimensional Viscoelastic Diffusive System

We examine a model for non-Fickian "sorption overshoot" behavior in diffusive polymer-penetrant systems. The equations of motion proposed by Cohen and White [SIAM J. Appl. Math., 51 (1991), pp. 472–483] are solved for two-dimensional problems using matched asymptotic expansions. The phenomenon of shock formation predicted by the model is examined and contrasted with similar behavior in classical reaction-diffusion systems. Mass uptake curves produced by the model are examined and shown to compare favorably with experimental observations

Dressed States of a two component Bose-Einstein Condensate

A condensate with two internal states coupled by external electromagnetic radiation, is described by coupled Gross Pitaevskii equations, whose eigenstates are analogous to the dressed states of quantum optics. We solve for these eigenstates numerically in the case of one spatial dimension, and explore their properties as a function of system parameters. In contrast to the quantum optical case, the condensate dressed states exhibit spatial behaviour which depends on the system parameters, and can be manipulated by changing the cw external field.Comment: 6 pages, including 6 figures. This paper was presented at ACOLS98, and is submitted to a special issue of J. Opt.

Algebraic Properties of Valued Constraint Satisfaction Problem

The paper presents an algebraic framework for optimization problems expressible as Valued Constraint Satisfaction Problems. Our results generalize the algebraic framework for the decision version (CSPs) provided by Bulatov et al. [SICOMP 2005]. We introduce the notions of weighted algebras and varieties and use the Galois connection due to Cohen et al. [SICOMP 2013] to link VCSP languages to weighted algebras. We show that the difficulty of VCSP depends only on the weighted variety generated by the associated weighted algebra. Paralleling the results for CSPs we exhibit a reduction to cores and rigid cores which allows us to focus on idempotent weighted varieties. Further, we propose an analogue of the Algebraic CSP Dichotomy Conjecture; prove the hardness direction and verify that it agrees with known results for VCSPs on two-element sets [Cohen et al. 2006], finite-valued VCSPs [Thapper and Zivny 2013] and conservative VCSPs [Kolmogorov and Zivny 2013].Comment: arXiv admin note: text overlap with arXiv:1207.6692 by other author

Scale-Free Networks are Ultrasmall

We study the diameter, or the mean distance between sites, in a scale-free network, having N sites and degree distribution p(k) ~ k^-a, i.e. the probability of having k links outgoing from a site. In contrast to the diameter of regular random networks or small world networks which is known to be d ~ lnN, we show, using analytical arguments, that scale free networks with 2<a<3 have a much smaller diameter, behaving as d ~ lnlnN. For a=3, our analysis yields d ~ lnN/lnlnN, as obtained by Bollobas and Riordan, while for a>3, d ~ lnN. We also show that, for any a>2, one can construct a deterministic scale free network with d ~ lnlnN, and this construction yields the lowest possible diameter.Comment: Latex, 4 pages, 2 eps figures, small corrections, added explanation

Development of a phase-change thermal storage system using modified anhydrous sodium hydroxide for solar electric power generation

A thermal storage system for use in solar power electricity generation was investigated analytically and experimentally. The thermal storage medium is principally anhydrous NaOH with 8% NaNO3 and 0.2% MnO2. Heat is charged into storage at 584 K and discharged from storage at 582 K by Therminol-66. Physical and thermophysical properties of the storage medium were measured. A mathematical simulation and computer program describing the operation of the system were developed. A 1/10 scale model of a system capable of storing and delivering 3.1 x 10 to the 6th power kJ of heat was designed, built, and tested. Tests included steady state charging, discharging, idling, and charge-discharge conditions simulating a solar daily cycle. Experimental data and computer-predicted results are correlated. A reference design including cost estimates of the full-size system was developed

Meanfield treatment of Bragg scattering from a Bose-Einstein condensate

A unified semiclassical treatment of Bragg scattering from Bose-Einstein condensates is presented. The formalism is based on the Gross-Pitaevskii equation driven by classical light fields far detuned from atomic resonance. An approximate analytic solution is obtained and provides quantitative understanding of the atomic momentum state oscillations, as well as a simple expression for the momentum linewidth of the scattering process. The validity regime of the analytic solution is derived, and tested by three dimensional cylindrically symmetric numerical simulations.Comment: 21 pages, 10 figures. Minor changes made to documen

Uniqueness of the electrostatic solution in Schwarzschild space

In this Brief Report we give the proof that the solution of any static test charge distribution in Schwarzschild space is unique. In order to give the proof we derive the first Green's identity written with p-forms on (pseudo) Riemannian manifolds. Moreover, the proof of uniqueness can be shown for either any purely electric or purely magnetic field configuration. The spacetime geometry is not crucial for the proof.Comment: 3 pages, no figures, uses revtex4 style file